Measuring the deviation from maximal mixing of 
atmospheric muon neutrinos at INO 
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• Neutrino factfile and D = 1/2 — sin 2 #23 

• u^, survival probabilities and binned up-down asymmetries 

• Error analysis and sensitivity to D 

This talk is based on my work [1] with S. Choubey. 

• Neutrino factfile 

We take a CPT-conserving framework with only three active light neutrinos: v e ,^,T 
as flavor eigenstates and ^1,2,3 as mass eigenstates with real, positive mass eigenvalues 



mi 2 3 and use A 



ij 
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m\. We relate the two sets of eigenstates by 



Here Upmns is defined in the standard CKM parametrization with possible Majorana 
phases ignored. We define 

D = 1/2 — sin 2 #23- (2) 
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mi\ < 0.057 eV 
30° < #12 < 38° 
36° < 23 < 54°, i.e. -0.16 < D < 0.16 
#13 < 12° 
5 = ? 



Table 1. Known 3a limits on neutrino parameters 

The known 3a limits on the relevant neutrino mass and mixing parameters are given 
in Table 1. It is also instructive to keep in mind that —0.10 < D < 0.10 at 2a and 
-0.07 < D < 0.07 at la levels. 



• Vp, survival probabilities and binned up-down asymmetries 

We explore the measurability of D from atmospheric muon neutrino and antineu- 
trino studies in a large magnetized iron calorimeter like ICAL - as proposed [2] in 
the India-based Neutrino Observatory INO. Effects due to matter and sub leading os- 
cillations on account of the solar neutrino mass scale A21 cause departures from the 
simple sin 2 2^23 dependence of 2-flavor vacuum oscillations in the survival probability 
Pfifi (Pp.fi) for a (Pjj). Both these effects as well as vacuum oscillation terms play 
a role in the measurement of \D\. However, matter effects are found to be crucial in 
resolving the octant ambiguity sin 2 #23 1 — sin 2 ^ 2 3 ; i-e. whether #23 is greater or 
less than tt/4, through a determination of the sign of D. Our numerical code exactly 
solves the three generation neutrino equations of motion in earth matter with the 
PREM density profile. We assume normal neutrino mass ordering, but we can do the 
analysis for the inverted case. We fix some of the oscillation parameters at chosen 
benchmark values, to wit A 31 = 2 x 1(T 3 eV 2 , A 2i = 8 x 1(T 5 eV 2 , sin 2 12 = 0.28, 
sin 2 2813 = 0.1, 5 = 0. Using the atmospheric neutrino flux of Honda et al [3] and the 
DIS cross sections of the CTEQ collaboration [4], we have simulated 14000 (7000) 
(Vfj) events for 1 MtonY of exposure for ICAL assuming a 50% trigger efficiency. We 
have binned our data both in energy E and the cosine of the zenith angle £. The dif- 
ference between the and up-down asymmetry ratios Un/D^ — Ua/ 'Da is shown 
in Fig. 1 with solid (dashed) pink and black lines corresponding to sin 2 #23 = 0.36 and 
0.5 respectively in matter (vacuum). The largest usable effects are in the topmost 
right panel (the lowermost right panel having problems due to earth-core uncertain- 
ties) with the 5-7 GeV energy bin corresponding to a maximum SPMAX in P MM . 
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Fig. 1. The difference U^/ D N — Ua/P>a for various bins of E and |c^|. 



• Error analysis and sensitivity to D 

We have done a \ 2 analysis using the pull approach, as detailed in [1]. Fig. 2 
shows the regions of sin 2 # 13 for which a maximal #23 can be rejected at 3cr (yellow), 
2cr (green) and la (magenta) levels. The capability to resolve the octant ambiguity 
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Fig. 2. Sensitivity to deviation from maximal 823 in terms of #13 

in 623 is shown in the A% 2 vs. sin 2 623 plot of Fig. 3 where A% 2 is the difference 
in x 2 f° r 'true octant' and 'false octant' values. The latter is made possible by the 
unique capability of a large magnetized iron calorimeter to isolate matter effects by 
distinguishing between and events. 
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Fig. 3. Capability of octant ambiguity resolution 
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